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1 .0  INTRODUCTION 


A  model  of  Iceberg  deterioration  must  attempt  to  use  the  best  available 
formulations  of  Ice  melt  In  sea  water  while  allowing  for  a  high  degree  of 
uncertainty  In  nearly  all  other  parameters.  This  uncertainty  starts  with 
the  size  and  shape  of  the  Iceberg  Itself.  The  Iceberg's  size  and  shape  are 
reported  by  an  observer  In  a  aircraft  flying  at  various  altitudes,  under 
varying  conditions  of  visibility,  and  at  widely  differing  ranges  to  the 
Iceberg.  Only  one  side  of  the  Iceberg  is  typically  seen.  Length  and  height 
estimates  are  approximate.  Nothing  of  the  underwater  shape  or  form  is  known 
to  the  casual  observer,  although  a  number  of  sonar  studies  have  been  conducted 
of  the  underwater  portion  (Robe,  1975,  and  Russell  et  al.,  1977).  These 
reports  detail  the  underwater  characteristics  of  a  small  sample  of  Icebergs 
and  from  this  sample  Idealized  Icebergs  can  be  constructed. 

To  this  uncertainty  of  the  shape  and  size  must  be  added  the  uncertainty 
of  the  determination  of  environmental  forces  which  account  for  the  deteriora¬ 
tion.  In  this  model  we  plan  to  use  only  three  envl ronmental  factors;  sea 
surface  temperature,  wave  height  and  period,  and  wind  speed.  The  almost  total 
lack  of  reporting  weather  stations  In  the  Labrador  Sea  and  on  the  Grand  Banks 
results  In  the  data  on  these  envl ronmental  parameters  being  extracted  from 
large-scale  computer  simulations  which  are  accurate  In  the  large-scale,  but 
lack  the  detail  often  needed  for  deterioration  computations. 

The  model  developed  here  does  not  attempt  to  be  more  rigorous  than 
the  data  justifies  or  that  requirements  demand.  It  draws  freely  on  rules 
of  thumb,  generalities,  approximations,  and  experience;  without  apology. 

From  the  uncertain  data  Input  for  the  Iceberg  characteristics  and  environ¬ 
mental  parameters  a  judgement  must  be  made  as  to  what  the  Initial  mass  of 
the  Iceberg  should  be  at  the  initial  point  of  the  deterioration  model. 

The  approach  In  this  report  Is  more  pragmatic  than  that  used  by  Kollmeyer 
(1966)  although  the  use  of  turbulent  heat  transfer  is  similar.  Kollmeyer 
constructed  a  heat  budget  for  the  Iceberg  which  balanced  the  heat  need  for 
melting  with  heat  sources. 


2.0  ICEBERG  CLASSES 


Icebergs  have  roughly  been  separated  into  types  or  classes  based  on 
form.  These  classes  as  presented  by  Murray  (1968)  are  given  In  Table  1. 

Table  I 

Iceberg  Shape  Classes 

(a)  BLOCKY  -  Steep  precipitous  sides  with  horizontal  or  flat  top. 

Very  solid  berg.  Length-height  ratio  less  than  5:1 

(b)  DRYDOCK  -  Eroded  such  that  a  large  U-shaped  slot  Is  formed,  with 
twin  columns  or  pinnacles.  Slot  extends  into  the  waterline  or 
close  to  It. 

(c)  DOME  -  Large  smooth  rounded  top.  Solid  type  berg. 

(d)  PINNACLED  -  Large  central  spire  or  pyramid  of  one  or  more  spires 
dominating  shape.  Less  massive  than  dome-shaped  berg  of  similar 
dimensions. 

(e)  TABULAR  -  Horizontal  or  flat-topped  berg  with  length -height  ratio 
greater  than  5:1. 

(f)  BERGY  BIT  -  A  mass  of  glacial  ice  smaller  than  a  berg  but  larger 

than  a  growler,  about  the  size  of  a  small  cottage.  Small  berg  or 

large  growler  Is  preferred  usage. 

(g)  GROWLER  -  A  mass  of  glacial  Ice  that  has  calved  from  a  berg  or  is 

the  remains  of  a  berg.  A  growler  has  a  height  of  less  than  three 

meters  and  a  length  of  less  than  six  meters. 

These  categories  are  arbitrary  since  In  reality  the  shape  classes  strongly 
overlap  each  other  (Farmer  and  Robe,  1977).  They  do  however  describe  certain 
gross  characteristics  of  Icebergs  which  relate  mainly  to  stability.  In  this 
report  we  limit  our  consideration  to  four  types  which  combine  the  character¬ 
istics  of  the  above  listed  classes.  The  four  types  are  defined  using  the 
classes  in  Table  1:  (1)  the  tabular,  which  combines  the  tabular  and  blocky 

classes;  (2)  the  domed,  which  combines  the  domed  and  the  growler;  (3)  the 
pinnacle;  and  (A)  the  drydock  which  are  the  same  as  In  Table  I. 
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3.0  DETERIORATION  OF  UNDERWATER  AREA  FOR  ICEBERG  CLASSES 


In  order  to  melt  the  underwater  surface  of  an  Iceberg  by  turbulent 
heat  transfer  one  needs  to  know  not  only  the  rate  at  which  heat  passes 
across  the  Ice/water  Interface,  but  also  the  surface  area  across  which 
the  heat  Is  transferred.  Since  an  observer  sees  only  the  shape,  height, 
and  width  of  the  above  water  portion,  the  task  of  deriving  the  underwater 
surface  area  can  neither  be  straightforward  nor  exact.  Since  we  are  trying 
to  derive  surface  area  I  have  chosen  to  relate  the  underwater  surface  area 
to  the  square  of  the  mean  waterline  length  so  that  the  result  Is  a  dimen¬ 
sionless  ratio. 

The  Iceberg  types  will  be  treated  as  stereotypical  constructs  from 
which  deterioration  factors  can  be  extracted  using  above  water  dimensions 
alone.  The  deterioration  factors  which  are  of  importance  are  waterline 
length  exposed  to  wave  erosion  and  the  underwater  surface  area  which  is 
exposed  to  forced  convective  melting.  The  waterline  length  Is  obtained 
directly,  but  the  underwater  surface  area  is  derived  from  the  length, 
width,  and  height  of  the  above  water  portion  of  the  iceberg  and  depends 
on  the  assumptions  which  are  defined  In  the  following  discussion  of  the 
stereotypical  Iceberg  types. 

a.  Tabular  Iceberg  Type 

In  Figure  1  a  perfectly  regular  rectangular  Iceberg  Is  portrayed 
with  a,  b,  and  c  being  the  width,  length,  and  thickness,  respectively.  The 
above  water  height  Is  represented  by  h  which  Is  equal  to  1/8  c.  For  a  blocky 
Iceberg  we  let  a  ■  b. 
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FIGURE  1.  An  Idealized  representation  of  an  iceberg  of  the  tabular  type. 

The  length  is  represented  by  b,  the  width  by  a,  and  the  total 
thickness  by  c.  The  above  water  height  h  is  equal  to  one  eighth 
of  the  thickness. 

The  hydrostatic  equation  requires  that  one  eighth  of  the  iceberg  be  above 
water  level  in  salt  water.  Therefore,  the  below  water  surface  area  A$  of 
this  iceberg  is  given  by 

As  =  ab  +  1.75ac  + 1.75  be  ( 1 ) 

The  average  waterline  length  squared  is  given  by 

Lw2  =  -5a2  *  25b2  +  .5as/a2  +  b2  < 2 ) 

Dividing  equation  (1)  by  equation  (2)  we  arrive  at  a  ratio  of  underwater 
surface  area  to  the  square  of  the  waterline  length  of 

Ratio: 


Let  b  *  Na,  then  (3)  becomes  ^ 

Ratio  = - 4N _ +  7(N+1)  /c\ 

2 * N?  -  2,/l  rNr  2  ♦  N2+  zj  1  *  N2  \) 
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The  ratio  of  thickness  to  width  c/a  normally  falls  In  the  range 

(Robe,  1975) 


Choosing  a  middle  value  of 


c/<s  =  .7 


(4)  becomes 


Ratio  t - 4N;  ♦  - (5) 

2  ♦  N*  ♦  2j  1  ♦  N2  2  *N2*  2/l  *  N2 

Therefore  the  underwater  surface  area  can  be  recovered  by  knowing  the  ratio  N 
of  the  longest  and  shortest  side  of  the  rectangular  iceberg  and  by  multiplying 
equation  (5)  by  the  square  of  the  observed  waterline  length.  This  result  will 
be  used  later  in  the  development  of  forced  convective  melting. 

b.  Pinnacle  Iceberg  Type 

We  now  follow  the  treatment  of  the  tabular  Iceberg  with  the  more 
complicated  double  pyramid  shape  in  Figure  2.  In  Figure  2  refers  to  the 
horizontal  dimension  at  the  largest  Iceberg  horizontal  section,  i'  Is  the 
waterline  length,  h  Is  the  vertical  half  thickness  of  the  Iceberg,  while  h' 

Is  the  above  water  height.  We  first  need  to  know  where  the  waterline  lies  on 
this  shape.  What  Is  the  value  of  h ' ?  Starting  with  the  fact  that  1/8th  of 
the  volume  Is  above  water  we  have  the  total  volume  V 

V  =  2/3/2h  (6) 

The  subaerial  volume  V  is 

a 

Va=V3n2^2nh  (7) 

where  and  h's  nh 
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FIGURE  2.  An  Idealized  representation  of  an  Iceberg  of  the  pinnacle  type. 

The  base  of  the  double  pyramid  Is  square  with  side  length  £  . 
The  altitude  of  each  of  the  double  pyramids  is  h.  The  above 
water  dimensions  are  height  h1  and  waterline  length 


The  above  water  volume  based  on  V8th  of  the  total  volume  from  (6)  Is 


VQ=  Vl2-£2h 


(8) 


Equating  (7)  and  (8) 

Vz  t?  s  ’/la 

or 


rt  =  .63 
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The  total  surface  area  of  the  double  pyramid  is  given  by 


A  =  4^2/4  +  h2 

and  the  submerged  area  Is  given  by 
As=3.2  +  h2 

The  average  waterline  length  Is 

Lw  =  '££lL±jL  =  1.21/  or 


.7Gj£ 


2 

The  squared  average  water  line  length  Is  then 

Lw2  =  58^ 2 


(9) 


(10) 


(11) 


The  ratio  of  the  submerged  surface  area  to  the  square  of  the  average  waterline 
length  is  given  by 

Ratio  =  -•52v/^2/4*h2 

/ 

Converting  to  £  and  |/ 

Ratio  =  5.52/25  +(hV)2 

A  reasonable  value  of  h//'  Is 
Is  chosen  the  ratio  of  the  submerged  area  to  the  square  of  the  waterline 
length  becomes 

Ratio  =3.2  (13) 


(12) 


(12)  becomes 


0.3,  Budd  (1977).  If  a  value  of  -  0.3 


The  submerged  surface  area  is  recovered  by  multiplying  the  ratio  by  the  square 
of  the  observed  waterline  length. 
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c.  Drydock  Iceberg  Class 

The  drydock  iceberg  is  the  most  stable  of  the  deteriorated  iceberg 
forms.  Time  is  required  to  cut  the  notch  through  the  center  portion  of  the 
iceberg  and  during  this  time  the  Iceberg  must  not  rotate.  The  drydock  berg 

Is  Idealized  in  Figure  3  where  the  width  Is  a,  the  thickness  c,  and  the  above 
water  height  b. 


FIGURE  3.  An  Idealized  representation  of  an  Iceberg  of  the  drydock  tvoe. 

The  ^eberg  Is  assumed  square  with  side  a  and  total  thickness  < 

IiJhS!a£$  *h1ch  for™Jthe  drydock  sides  and  the  drydock  cut  are 
each  taken  as  one  third  a.  The  height  of  the  above  water 
portion  Is  represented  by  b. 


A  berg  which  Is  drydocked  will  be  assumed  to  have  a  length  nearly 
equal  to  Its  width.  Following  the  arguments  In  the  preceding  sections  we  have 
from  Figure  3  that.  If  b  is  the  depth  of  the  drydock  cut,  the  above  water 
volume  Is 


Va=%a2b  (14) 

and  the  below  water  volume  Is 

V#=o2(c-b)  (15) 

If,  as  Is  required,  the  above  water  Ice  volume  Is  one  eighth  of  the  total 
volume  then  using  equations  (14)  and  (15) 

%a2b2  =  V7a2(c-b)  (16) 

or 

5.66b  *  c  (17) 

The  below  water  surface  area  then  becomes 


A$  =  a2  +  4o(c-b) 
and  substituting  from  (17) 

:  o2  ♦  4o  (4.66b) 

The  average  waterline  length  for  this  berg  Is 

Lw=  g±^?q—  or  i.21a 
2 

The  average  waterline  length  squared  becomes 
Lw2  =  1.46a2 


(18) 


(19) 


(20) 
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ratio  of  the  underwater  surface  area  to  the  square  of  the  average 
waterline  length  Is  given  by 

Ratio  =  .68  +  12.79%  (21) 

The  value  of  ^  can  easily  be  obtained  from  observation. 

If,  for  example  %=  .5  then 

Ratio  =  7.08  {u) 

The  underwater  surface  area  subject  to  forced  convection  can  then  be  found  by 
multiplying  this  ratio  In  (22)  by  the  square  of  the  observed  waterline  length. 

d.  Domed  Iceberg  Class 

This  Is  the  class  of  Icebergs  which  Is  more  or  less  spherical.  They 
are  so  because  they  are  In  an  advanced  state  of  deterioration  and  roll  over 
at  the  slightest  disturbance.  The  domed  class  Includes  the  smaller  bergs 
and  growlers  as  well  as  an  occasional  large  Iceberg  which  happens  to  have 
gone  spherical  early.  To  represent  the  domed  iceberg  a  sphere  was  chosen. 
Figure  4,  and  the  computations  were  made  using  spherical  coordinates.  The 
radius  Is  represented  by  p  ,  the  vertical  rotation  by  *  and  the  horizontal 
rotation  by  9  .  The  sphere's  center  lies  a  distance  d  below  the  water 
surface.  From  Figure  4  we  can  compute  the  above  water  surface  area  from 

% 

A  az  Jr2 sinddddd 

•  • 


Aa=  27rp2(1-<%) 


(23) 


The  Iceberg  is  assumed  spherical  with  the  sphere  having  a 
radius.  The  waterline  Is  a  distance  d  above  the  sphere 
center.  Spherical  coordinates  are  used  for  calculations. 


Subtracting  (23)  from  the  total  surface  area  of  a  sphere  one  obtains  the  below 
water  surface  area. 

A« *  4jrA2-2jrp2(l-d/p)  (24) 

Since  one  eighth  of  the  total  volume  Is  above  water  It  is  possible  to 
compute  d  or  the  distance  from  the  sphere  center  to  the  waterline  (Figure  5). 

The  shaded  portion  of  the  circle  Is  proportional  to  the  volume  of  that  portion 
of  a  sphere  represented  by  the  shaded  area.  The  area  Is  given  by 

Aa*2/jriy  (25) 

where  xs  Vi  ♦(y+d)2 
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so  that  (25)  becomes 


i-e _ _  , 

*a*2f  J-\  -(y>drdy 

let  u  *  y+d  and  du  ■  dy 
substitute  into  (26) 


integrate  and  substitute  back  to  get 

Aa  =  frk  -  d/2yi-d2  -  V2  Asintd) 


Equation  of  a  unit  circle 
(y+d)J  *  X1  :1 


FIGURE  5.  The  above  water  fraction  of  a  domed  Iceberg  Is 

calculated  by  integrating  the  volume  elements  xdy. 


The  area  of  a  circle  of  unit  radius  is  equal  to  t  .  The  1 /8th  portion  would 
then  be  n/&»  and  from  (28) 


7q  z  -  d>/Tdr  -  Asin(d)  (29) 

or  d?. 6347  (30) 

for  a  unit  sphere. 

Waterline  length  as  a  function  of  the  radius  is  then 

P2-(.6347P)2  :  (lw/2)2  (3D 

or 

Lw=  1.54P  (32) 

The  waterline  length  squared  becomes 

Lw2--2.39P2  (33) 

To  compute  the  ratio  of  the  submerged  surface  area  to  the  square  of  the 
waterline  length  we  use  (24)  and  (33). 

Ratio  :|4ffp2-2n-p2  (l-*%)J/2.39  (34) 

or  using  (30) 

Ratio  =  4.29  ,35, 

To  recover  the  underwater  surface  area  multiply  (35)  by  the  square  of  the 
observed  waterline  length. 
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4.0  ICEBERG  MASS  ESTIMATE 


The  total  mass  of  an  Iceberg  Is  difficult  to  obtain  except  by  using 
stereo  photography  and  a  great  deal  of  computational  work.  An  estimate 
of  the  mass  Is  necessary,  however,  as  an  Initial  point  for  a  deterioration 
scheme.  Farmer  and  Robe  (1977)  were  able  to  develop  an  extremely  simple 
method  for  estimating  the  mass  of  Icebergs  of  random  shape  and  size.  Their 
method,  which  Is  based  on  a  number  of  icebergs  measured  carefully  using 
stereo  aerial  photography,  depends  only  on  the  length,  width,  and  height 
of  an  Iceberg.  Their  relationship  gives  the  Iceberg  mass  M  In  metric 

3 

tons  (10  kg)  as: 

M  »  3.01  x  Length  x  Width  x  Height  (36) 

with  the  units  of  length,  width  and  height  all  In  meters. 

5.0  ENVIRONMENTAL  FACTORS 

Only  three  environmental  factors  are  considered  In  this  model:  (1)  wind, 
which  accounts  for  a  portion  of  the  berg  motion  relative  to  the  water;  (2)  sea 
surface  temperature,  which  will  provide  the  thermal  driving  force  for  melting; 
and  (3)  wave  height  and  period,  which  creates  a  very  turbulent  near  surface 
zone  eroding  the  Ice  rapidly  near  the  water  line. 

a.  Wind 

Wind  forces  applied  to  the  above  water  portion  of  the  Iceberg 
will  cause  the  Iceberg  to  move  relative  to  the  water  thereby  Inducing  forced 
convective  melting  of  the  Iceberg.  Icebergs  are  generally  In  motion  relative 
to  the  water  at  any  particular  depth,  even  In  a  zero  wind  situation,  because 
of  the  change  of  current  velocity  with  depth  and  the  fact  that  the  iceberg 
moves  In  response  to  the  Integrated  force  applied  to  It  by  water  at  all 
depths.  A  good  value  for  this  background  relative  motion  aV  between 
Iceberg  and  water  Is  0.05  m/s.  The  relative  motion  caused  by  the  wind 
was  given  by  White  et  al.  (1980)  as: 
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dVzW/O  +  B). 


(37) 


where  B  Is  shape  factor  for  different  icebergs  and  W  is  the  steady  wind 
speed.  The  shape  factor  was  estimated  as  follows: 

Pinnacle  berg:  B  *  44+5 
Dome  or  drydock:  B  *  66+8 
Tabular  berg:  B  »  75+9 


If  the  a  mid  value  of  B  *  65  is  chosen  for  equation  (37)  and  the  0.05  m/s 
background  relative  motion  is  added  in  (37)  becomes: 


4V  = 


W  +  3.30 
66 


(38) 


where  4V  and  W  are  in  m/s. 

b.  Sea  Surface  Temperature 

Sea  water  temperatures  in  the  region  around  the  Grand  Banks  of 
Newfoundland  are  highly  variable.  The  water  temperature  can  vary  many  degrees 
in  short  distances  and  over  a  few  meters  depth.  The  sea  surface  temperature 
can  either  be  higher  or  lower  than  the  surface  layer  which  surrounds  the  ice¬ 
berg.  However,  lacking  the  ability  to  sample  water  profiles  in  this  area,  sea 
surface  temperature  remains  the  best  available  Indicator  of  the  ability  of  the 
water  to  melt  an  iceberg.  The  difference  between  the  water  temperature  and 
the  melting  temperature  of  Ice  in  30°/oo  seawater  (approximately  -1.63  °C) 
provides  the  thermal  driving  force  for  melting  the  Iceberg  Ice. 

c.  Wave  Height  and  Period 

The  erosion  of  the  waterline  by  wave  action  Is  one  of  the  most 
forceful  mechanisms  of  Iceberg  deterioration.  The  deterioration  takes  place 
directly  by  wasting  and  indirectly  by  calving  of  ice  from  the  undercut  sides 
of  the  Iceberg.  Both  the  wave  height  and  wave  period,  together  with  the  water 
temperature,  play  a  large  role  In  developing  the  wave  cut  notch.  Short  period 
waves  tend  to  deepen  the  wave  cut  notch  faster  than  waves  of  longer  periods. 
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6.0  DETERIORATION  RATES 


The  turbulent  heat  transfer  between  the  seawater  and  the  Iceberg  will  now 
be  calculated.  The  transfer  of  heat  is  dependent  on  the  surface  area  exposed, 
the  seawater  temperature  and  the  relative  ice /water  speed.  The  loss  of  mass 
through  wave  undercutting  and  calving  will  also  be  considered. 


a.  Side  wall  meltinc 


Modifying  the  development  in  White,  et  al.  (1980),  the  rate  of 

melting  V  in  cm  s'1  °C"'  is. 

*  m 


where  9:  is  the  density  of  ice  (0.9  gm  cm  )t  r  is  the  latent  heat  of 
fusion  of  water  (344  J  gm"  ),  and  #wis  the  heat  transfer  rate  per  unit 

0  1  ft  1 

area  per  degree  C  in  J  cm"*  s"  °C~  * 


The  local  Nusselt  number  which  relates  the  local  heat  transfer  rate  to 
the  thermal  conductivity  is  given  by: 


where  x  is  the  downstream  dimension  in  (cm),7vvis  9^ven  above,  K  is  the 
kinematic  thermal  conductivity  Js~^  cm"^  °C"^  (we  will  use 
K*  .006  Js"^cm'^°C~*  which  Is  the  value  for  fresh  water  at  17.5  °C 
and  varied  only  slightly  with  the  temperature  and  salinity),  and  AT  1 s  the 
thermal  driving  force  or  water  temperature,  Tw,  minus  the  Ice  melting 
temperature,  T^. 

Presenting  the  Nusselt  number  In  a  different  form  where  it  is  based  on 
the  waterline  length. 


Nu^:  0.056Re^°-8Pr°-4 


(41) 
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Where  Ltt  Is  the  waterline  length  In  (cm),Re..1s  the  waterline  Reynolds 
number  ,(Re:  L°/ f AV  is  the  relative  iceberg/water  velocity 

(cm  s’  1),  is  the  kinematic  viscosity  of  sea  water  (cm^  s’* )  which 
for  water  of  30  °/oo  and  density  (Fw*  1.0275  gm  cm’3)  is  represented  by 


*w  =  0183-6.18x10  T  + 1.51x10  t2- 2.0x10  t3  (42) 

with  only  a  slight  salinity  dependence,  and  Pr  is  the  Prandtl  number  (**r:  / *J 
where  *  is  the  thermal  diffusivity.  The  Prandtl  number  can  be  represented  by 

Fr  - 1.13  -  0.431  T  -  0.000893T2*  0.00521  T3  (43) 


Combining  (40)  and  (41)  and  solving  for®  we  obtain: 


9w  ■-  0.056-^- 

Lo 


substituting  in  (39) 


.  0.056  /  K4T\/AVlo\08  o.4 

[  Lo  Aw 


>.006x0,056 

0.9x344 


if  the  freezing  point  for  ice  in  30° /oo  seawater*! is  equal  to  -1.63°C  then 
(converting  centimeters  to  meters) 


vm  »  ^09^0a(\;L63)(^-)pr°< 


in  ms'1  V1 


When  computing  AV  use  the  wind-driven  relative  velocity  plus  0.05  m/s  which 
accounts  for  the  relative  velocity  of  the  Iceberg  and  water  In  a  zero  wind 
situation. 
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b.  Cal vl ng 


The  direct  melting  by  wave  turbulent  heat  transfer  at  the  waterline 
will  not  be  treated  separately,  but  will  be  lumped  In  with  the  problem  of 
calving  as  a  result  of  wave  undercutting.  Calving  occurs  when  an  undercut 
slab  of  Ice  on  the  sidewall  of  an  Iceberg  breaks  off  under  Its  own  weight. 

The  frequency  with  which  this  happens  depends  upon  wave  height  H,  wave  period 
t,  and  the  thermal  driving  force.  The  amount  of  Ice  that  Is  lost  depends  upon 
the  shape  of  the  iceberg,  the  length  of  the  face  exposed  to  waves  and  to  the 
height  of  the  calving  face.  An  assumption  Is  made  that  the  roughness  of  the 
iceberg  Is  nearly  1%  of  the  wave  height,  i.e.  2  cm  surface  roughness  for  waves 
in  the  1  to  3  meter  range.  From  White  et  al .  (1980)  the  speed  with  which  the 
water  line  is  cut  Is  then 

Vw  *  0.000058 H/t  (47) 

with  Vw  being  the  erosion  velocity  in  ms"1  °C"1 ,  H  the  wave  height  in 
meters,  and  t  the  wave  period  in  seconds. 

Calving  which  results  from  this  undercutting  Is  assumed  to  take 
place  along  a  front  which  is  one  half  of  the  observed  waterline  length. 

This  factor  Is  chosen  because  the  wave  activity  tends  to  concentrate  at 
a  few  places  along  the  waterline  and  In  any  case  only  one  side  of  the 
Iceberg  at  a  time  is  exposed  to  the  wave  erosion.  Rarely  does  the  entire 
side  of  an  iceberg  calve  at  one  time.  The  amount  of  Ice  calved  depends  not 
only  on  the  length  of  the  calving  wall  but  also  on  the  wall  height  or  the 
calving  slab  thickness.  For  a  given  slab  thickness  above  a  wave  cut  notch 
a  criterion  has  to  be  developed  for  when  the  notch  has  been  cut  deeply  enough 
Into  the  side  of  the  Iceberg  to  cause  calving.  However,  first  a  guideline  for 
the  calving  slab  thickness  must  be  stated  as  a  function  of  the  total  above 
water  height.  For  tabular  and  drydocked  icebergs  with  their  fairly  vertical 
sides  a  value  of  80%  of  the  observed  height  Is  chosen  from  observation.  For 
pinnacle  Icebergs  with  their  more  sloping  side  wall  a  value  of  30%  of  the 
Iceberg  above  water  height  Is  chosen.  Domed  Icebergs,  because  their 
instability  leads  to  frequent  changes  In  waterlines,  are  not  permitted 
to  calve  In  this  model. 
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The  critical  wave  cut  waterline  notch  depth  which  gives  rise 
to  calving  comes  from  White  et  al .  (1980)  and  is  given  by: 

JLT  '•  0.33(3.75H  ♦  Sj2)1/a  (48) 

where  H  is  the  wave  height  in  meters  and  (S^)  is  the  slab  thickness  in 
meters.  When  the  depth  of  the  wave  cut  notch  developing  at  a  velocity  Vw 
given  in  (47)  equals  or  exceeds  £r  in  equation  (48)  calving  occurs.  In 
addition  the  dimensions  of  the  iceberg  are  reduced  by  the  value  £r. 

7.0  THE  MODEL 

The  elements  of  a  practical  iceberg  deterioration  model  have  been 
constructed  in  the  proceeding  sections  of  this  report.  The  application 
then  becomes  a  matter  of  using  the  formula  developed  to  estimate  the  speed 
with  which  the  Iceberg  Is  destroyed.  The  steps  are  as  follows,  and  in 
Figure  6: 


a.  Collect  the  data  on  size  and  shape  of  the  iceberg  together  with 
the  wind  speed,  sea  surface  temperature,  and  wave  height  and  period  for 
the  iceberg's  location.  (Source:  International  Ice  Patrol,  Governors 
Island,  New  York,  and  FNOC,  Monterey,  California.) 

b.  Determine  the  initial  mass  of  the  Iceberg  in  metric  tons  using 
equation  (36). 

c.  For  the  appropriate  iceberg  class  compute  the  underwater  surface  area 
from  the  square  of  the  waterline  length  using  equations  (5),  (13),  (22),  and 
(35). 


d.  From  equation  (38)  compute  the  speed  of  the  Iceberg  relative  to  the 
water  using  the  wind  speed. 

e.  Enter  the  speed  of  the  Iceberg  relative  to  the  water,  the  waterline 
length,  and  the  sea  surface  temperature  Into  equation  (46)  to  give  the  melting 
speed  In  m  s“*  °C“*.  Multiply  this  result  by  the  thermal  driving  force, 

the  sea  surface  temperature  plus  1.63.  Multiply  this  In  turn  by  the 
underwater  surface  area.  This  provides  the  forced  convective  melting  in 
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ms".  Then  multiply  this  result  by  the  deterioration  time  available  in 
seconds  and  by  the  density  of  ice  to  obtain  the  mass  lost  in  metric  tons  by 
forced  convective  melting. 

f.  Use  the  time  available  for  deterioration  together  with  the  sea 
surface  temperature  and  wave  height  and  period  in  equation  (47)  to  determine 
the  extent  of  wave  undercutting.  If  this  value  exceeds  the  critical  value 
from  equation  (48)  allow  the  iceberg  to  calve  along  1/2  of  the  iceberg 
length.  If  the  undercutting  exceeds  twice  the  critical  value  in  equation 
(48)  allow  the  iceberg  to  calve  twice,  etc. 

g.  Subtract  the  mass  lost  by  forced  convective  melting  and  calving  from 
the  original  mass  and  shorten  the  iceberg  by  the  length  of  the  wave  cut  notch. 

h.  Compute  a  new  underwater  surface  area  and  repeat  the  calculations  for 
a  new  time  period. 
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*  DOMED  ICEBERG  TYPES  ARE  NOT  PERMITTED  TO  CALVE. 


FIGURE  6.  Flow  diagram  for  preliminary  Iceberg  deterioration  model 
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8.0  VERIFICATION  AND  TESTING 


A  model  of  a  physical  real-world  system  must  be  verified  and  tested  In 
order  to  be  practical.  Testing  of  the  model  must  account  for  the  totality  of 
the  deterioration  effects  and  not  each  detail.  The  Important  outcome  Is 
whether  the  Iceberg  as  a  whole  melts  on  time  and  not  whether  each  particular 
Input  to  the  deterioration  Is  exact.  The  verification  can  only  be  statistical 
for  what  in  essence  Is  a  statistical  process. 

Since  verification  and  testing  of  the  model  must  proceed  from  the 
full-scale  and  from  the  open  ocean  environment,  the  testing  Is  likely  to  be 
difficult.  A  statistical  sample  Is  required.  The  sample  size  has  to  be  large 
enough  to  be  significant,  but  small  enough  to  be  managable.  A  sample  of  the 
order  of  10  for  each  Iceberg  type  Is  considered  satisfactory. 

A  verification  experiment  should  track  ten  to  fifteen  Icebergs  of  each  of 
the  four  types  from  the  time  they  reach  open  water  outside  the  Ice  pack  until 
they  entirely  melt.  The  measurements  collected  for  the  individual  icebergs 
must  correspond  with  the  measurements  which  would  actually  be  used  in  the 
model .  The  measurements  shoul d  be  repeated  at  two  or  three-day  1 nterval s 
for  the  life  of  the  Iceberg. 

Measurements  are  critical  during  the  rapid  deterioration  that  takes  place 
In  water  warmer  than  5°C  and  during  stonpy  Intervals. 

As  an  Initial  phase  of  testing,  those  Icebergs  In  the  southern  portions 
of  the  Ice  Patrol  area  which  are  undergoing  the  most  rapid  melting  should  be 
followed. 
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